Hydrodynamics of an electron-positron plasma near a black hole; applications to jet 

formation. 
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We investigate some features of the hydrodynamics and neutrino physics in the (predominantly) 
electron-positron plasma above a hyperaccreting disk (or torus) around a black hole, a conjectured 
engine for a short gamma ray burst. We suggest a possible scenario in which plasma in the region 
very near the black hole, energetically driven by neutrino annihilation, emerges as a subsonic wind, 
which in a spherically symmetrical case would decelerate as it moves out. In this case we argue that 
the plasma heating will be primarily through neutrino-electron and neutrino-positron scattering, and 
that this process will be important throughout a region considerably larger than that of the neutrino 
annihilation process. In simple solutions a relatively gentle anisotropy in the heating through this 
process can create an approximately conical sonic surface, aligned with the system's axis. Inside 
this cone the fluid accelerates upwards as in standard jet models. 

PACS numbers: 98.70.Rz, 98.62Mw 



INTRODUCTION 

A popular model for the mechanism of short gamma 
ray bursts begins with a "hyperaccretion disk" (or torus) 
created in the last moments of a neutron star-black hole 
merger An electron pair plasma is then created 

in a region of relatively low density along the disk axis 
through the process v + v — > e + + e~, where the neutri- 
nos originate in the dense inner part of the torus. This 
pair production is concentrated in a region very near the 
black hole, since it is this region that combines proxim- 
ity to the neutrino production region with the geometric 
advantage that is given when neutrino beams originating 
from opposite sides of the torus meet head-on in the re- 
gion of the axis. Indeed, in the simulations reported in 
ref. [H the energy density deposited in the plasma per 
unit time near the axis is proportional to where z is 
the distance to the center of the black hole. 

The plasma is then to be accelerated, by hydrodynam- 
ics alone, to form aT~ 100 jet as it moves out the disk 
axis. 

In the case of high rates of accretion, the densities and 
temperatures of the matter in the dense torus render the 
matter opaque, or somewhat opaque, to neutrinos in the 
innermost region. Thus both the disk dynamics and the 
proposed production mechanism depend on a neutrino 
transport calculation. Depending on choices of input pa- 
rameters, as well as on the computational framework for 
estimating production rates, one may arrive at the con- 
clusion that the model is sufficient to supply the jet for 
an energetic gamma ray burst, or at the conclusion that 
it is insufficient. 

In the present note we propose an alternative to this 
picture, while adhering to the general description that an 
plasma formed on the periphery of the torus can be 
accelerated by hydrodynamic forces to form an extremely 
rclativistic jet along the direction of the disk axis. Taking 
a three solar mass black hole, as in the simulations of ref. 



Q, with a Schwarzchild radius of approximately 7 km., 
let us circumscribe the system with a sphere of radius, 
say, 50 km. and ask what happens in this region. In 
the models used in ref. [H , the background density at this 
distance is fairly small, even at polar angles 9 up to 60 
deg.; in model B of this reference it ranges from about 
10 6 g cm -3 along to axis to about 10 8 g cm -3 for 9 =60 
deg. 

We shall describe a scenario in which this whole re- 
gion fills with e plasma, with small outflow velocity, 
in which the variation of plasma temperature with po- 
lar angle drives a subsonic— > supersonic transition along 
a more or less conical surface, as shown in fig.l. Look- 
ing at fig.l we should state at the outset that some of 
the features are calculated later in this paper, and some 
imagined. Two key elements that are calculated or esti- 
mated are the exhibited ratios of temperatures and the 
shape of the sonic transition surface. 

The figure shows a slice of the system by a plane that 
is perpendicular to the plane of the accretion disk and 
which contains the axis of the disk. The cross-sections 
of an inner torus and that of an outer torus are shown 
with different shadings. In the inner one we expect the 
neutrino opacity to be dominated by scattering and ab- 
sorption on nucleons, as it is in the simulations of refs Q 
-0- Under the conditions of energy output that we con- 
sider, the innermost torus may be v opaque; it may be 
marginally transparent. The nucleon density decreases 
steadily as one moves from the center to the edge of this 
torus 1 In the surrounding toroidal region, less shaded, 



1 Most theories of the torus (or disk) proper, e.g. ref. [5j assume 
relatively constant baryonic density as one moves up from the 
median plane to the surface, whereas if one demands approximate 
hydrostatic equilibrium in the presence of the (tidal) force needed 
to sustain rotation around a point on the axis, but not at the 
center of the system, then the baryonic density must decrease 
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the considerations of ref. 11 1 




T=5 MeV 



T=8 MeV 




FIG. 1: Cross-section of the accretion torus, showing the black 
hole at the center. In the inner torus with darker shading, v- 
nucleon reactios dominate the v opacity. In the lighter gray 
area e + , e~ plasma dominates. The edge of the lightly shaded 
area is our injection sphere. The curved lines at the top form 
the boundary between subsonic and supersonic flow, and the 
arrows indicate flow of e + ,e~ plasma across this boundary. 
The temperatures are given in MeV. 



as we move from the inside to the outside, the leptonic 
opacity becomes significant. In the more lightly shaded 
torus leptonic scattering takes over completely. The very 
lightly shaded region is occupied by the plasma halo of 
our model. The v opacities are small enough so as not to 
change the v fluxes and energy distributions at the sur- 
face of the circumscribed sphere, our "injection" sphere. 

In estimates that we make in more detail in appendix 
A we find: 

a) That for a neutrino-surface temperature of 5 MeV, 
located approximately at the edge of the inner (dark) 
torus, the average temperature of the overlying plasma 
is maintained at a temperature of about 3 MeV, with 
the heating provided by neutrino-electron and neutrino- 
positron collisions, and the cooling being the usual e + + 
e~ — > v e + v e process. This heating mechanism was sug- 
gested by Woosley jlfj and also plays a central role in 



drastically as one approaches the boundary of the disk. This 
point was made in the last section of ref. It could be 

relevant to the model being proposed here, since it provides a 
picture of a smooth transition between the disk interior and the 
pair-plasma atmosphere being proposed here. 



b) That the v optical depth of our plasma sphere is 
approximately .1 under the above conditions. 

In slightly more realistic estimates with the tempera- 
ture varying with the distance from the core we find that 
the region outside the large sphere contributes only a 
small amount of additional optical depth; this was be- 
hind our choice of 50 km as the radius. All of this 
begs the question of how the plasma got there in the 
first place. We have made estimates based on seeding 
through the v e + D e — > e+ + e~ process, followed by 
heating through the scattering processes, that indicate 
a development time of the order of .01 sec. 

We believe that there is an internal consistency in our 
assumptions, but we cannot know that above picture will 
be sustained by future detailed calculations of the com- 
plete accretion torus cum transport calculation. We do 
believe that the simulations to date lack an essential piece 
of physics in the above speculations, namely the seeding 
of the far out region with some pair plasma followed by 
plasma build-up through scattering processes. 

Since our picture has at least an order of magnitude 
more energy deposited by the neutrinos that leave the 
torus core than does the conventional model, we have 
an incentive to look for acceleration mechanisms that 
can convert some appreciable fraction of the energy to 
a relativistic jet. The purpose of the present paper is to 
explore the hydrodynamics of the plasma configuration 
described above. Our main result, which admittedly will 
feed in more assumptions, will be that the temperature 
variation of the plasma as a function of 9, coming from 
the greater proximity of the torus for large 9, at a given 
radius, can at the same time enable the passage from 
subsonic to supersonic flow speeds and concentrate the 
flow in a cone along the axis. 

The calculation assumes a steady flow of plasma across 
the surface of the large sphere, at a subsonic speed, as 
the boundary condition for a steady flow solution of the 
equations of fluid dynamics. However, for orientation, 
in the next two sections we consider examples of several 
kinds of flows, both time-independent and time depen- 
dent. 



ISOTROPIC INJECTION IN THE REGION 
CLOSE TO THE HORIZON 

As we shall consider only the temperature domain T > 
2 MeV in which the electrons are quite relativistic, we 
take the pressure to be given by P = bT 4 /3, where b = 
(ll7r 2 /60) for an electron-positron-photon plasma, with 
negligible contribution from nucleons. For the energy 
density, in some later examples, we add the contribution 
of nucleons of number density n, p — bT 4 +nM, where M 
is the nucleon mass. We will remain in domains in which 
the contribution of nucleons to the pressure is negligible. 
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The perfect fluid equations are 



dP 



A1V 



1 d 



-Wg(P + p)U»U»] 



r» x u v u x (p + P ) = ft". 



In the case in which n ^ the equations ([T]) are sup- 
plemented by the equation of baryon conservation. Here 
ft" has non- vanishing components ft , ft 1 * that are the re- 
spective time rates of deposition of energy and momen- 
tum density in the plasma by neutrino annihilation, as 
seen by an observer at infinity. We comment more on 
these functions below. Wc define 



/0) = -ffooO) = (1 - r s/r) , 
and introduce the variable y of [l3[ and p2 
y=[f(r)/(l-v 2 )} 1 / 2 , 



(2) 



(3) 



Henceforth we choose the Schwarzschild radius as the 
unit of distance, r s — 1. Specializing to the case of n = 0, 
P = p/3, and substituting for the four- velocities in ([TJ) , 



f(r) 



-1/2 



vf(r) 1 / 2 



we obtain, after much algebra, 



a 



a 



(4) 



(5) 



and 



Or 



log[py 4 ] 



(1 - v 2 )py 2 
dp 



—v 



at 
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a 



vfh°) 



dt\l-v 2 



(6) 



The neutrino annihilation sources that we simulate 
with ft" in (|5|) and ([6|) produce pairs nearly at rest in 
the system of the observer at infinity, and at a rate that 
is essentially independent of the flow of the plasma into 
which they are injected. We implement this picture by 
adding an amount of energy density Ap, measured in an 
inertial frame, at rest, at some point in the vicinity of the 
hole, in time, At, as measured by the observer at infinity. 
Then the right hand side of the eq Q for A(<9 p T"°) is, 



ft = A(c^T"°) = [U°(v = 0)} 2 s = f- 



s, 



(7) 



where the four velocity U° is given by and s is the 
energy density deposition rate, s = (Ap/At). We take 
ft r = 0. 

Since other authors use different source terms for the 
same equations, we elaborate a little on our choice. If in 
the case ft" = we move everything in ([T]) except the time 
derivative terms to the right hand side, then we describe 



the equations as writing the change in time of T 00 and 
T 0r , respectively, as a sum of changes due to compression 
and due to transport. When we introduce the neutrinos, 
we add the additional changes to the energy-momentum 
tensor coming from their annihilation. This change, over 
the infinitesimal AT, is completely independent of the 
fluid motion. 

In contrast, in the case of the source terms use d by 
Thompson, et al. [HI], as well as those used by ref. [lq . 
and, we believe, by ref. [H, induce changes in T 00 and 
T 0r that depend on the velocity of the fluid. They are the 
correct equations when the source arises from radiative 
transfer, but not appropriate for our case in which AT° : ° 
and AT°' r = are to be calculated by the observer at 
infinity. 

It is nevertheless instructive to calculate the source 
function when we assume that energy density (and no 
momentum) is deposited at a given rate, q, in the co- 
moving system. The sources ft and ft 1 * are then, 



h r = 



h° = 
Ap 



Ap 



[U°(v)f = J 



At 



At. 

U°(v)U r (v) =vqy. 



(8) 



where we have now defined q to be the rate of energy 
density deposition as measured by an observer in a co- 
moving inertial frame, 



Ap 



AtVT^VJ 



(9) 



Note that in this case the combination that provides the 
source term in (|6|) vanishes, ft r — vfh° — 0. 

For the case of steady flow, the equations (HJ) are equiv- 
alent to those of Thompson et al [15| , the latter special- 
ized to the case of vanishing baryon density. 2 

We begin discussion of the solutions with the case in 
which we set the right-hand sides of both ([5]) and ^ 
equal to zero, looking for static solutions in the case of 
energy injection through a bounding surrface, but with 
no volume injection of energy. These equations now are 
exactly those of Flammang [3] (see also refs. [HI and ref. 
[I3|) , specialized to the case of the completely relativistic 
plasma with no nucleons present. Picking boundary con- 
ditions such that a constant temperature is maintained 
at an injection radius r$ > 1, and injecting at this point 
with velocity vq we find, 



= (wo - «o) 



r 2 f[ro) 



(10) 



2 The equations of ref. Il5l derive in turn from those of ref. |16|| . 
further specialized to include only the first moment term in the 
energy absorption part of the rediative transfer equation. 
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In the P = p/3 limit the magnitude of the energy density 
(i.e temperature) at the injection surface scales out of the 
velocity determination. 

Differentiating (fTT)|) we have, for all points above the 
injection radius, 



dv 
dr 



(l-3^ 2 ) 



dv 
dr 



(l-v^/v*) 



f(r) 



(11) 



where v s is the speed of sound. Noting that the right 
hand side of (fTTj) vanishes only at the point r = 3/2, we 
conclude that when the injection radius is outside this 
point the flow cannot make a transition from subsonic 
to supersonic (or vice- versa); it is easy to see that sub- 
sonic injection leads to a flow that decelerates, even in 
the absence of gravity, as it moves upward, gradually fill- 
ing all space with plasma and trading kinetic energy for 
internal energy. The velocity approaches zero at infinity. 
By contrast, supersonic injection into a surface outside 
r = 3/2 leads to an accelerating flow, with internal en- 
ergy being traded for kinetic energy and an asymptotic 
velocity of c. This is the situation envisioned in models 
of a gamma ray burst jet. (All of this, of course, is for the 
perfect P = p/3 plasma; in a real case, the plasma even- 
tually cools to the point at which the e^- are no longer 
relativistic.) 

The options change when the plasma is launched from 
below the radius r = 3/2. Now there can a transition 
from subsonic (below) to supersonic (above) for upmov- 
ing plasma if the injection velocity has exactly the right 
value (a function of ro) at the injection surface. If the 
injection velocity vq is greater than this critical value, 
there is no physical solution. If the vq is less than the 
critical value, then the velocity will peak at r = 3/2 and 
decrease thereafter. Fig. 2 shows plots of the solutions to 
([TO]) both for a substantially subcritical injection speed 
and for a case that is subcritical by a tiny amount, the 
latter plot showing the non-avoided crossing of two roots 
that is the key to the transition. 

All of the above analysis, and in particular the plots of 
fig. 2, are essentially the same as presented by Flammang 
14} . In the application to accretion flows in this reference 
and others, where v < 0, the transition is from subsonic 
(above) to supersonic (below), and the subcritical solu- 
tions represented by the (dashed) curves in fig. 2 can 
be dismissed as unphysical 3 . In our problem, though, 
we shall argue that the solution analogous to the lower 
dashed curve in fig. 1 could be the one chosen, once we 
put in the source terms, h° 7 h r . 



3 In the case of accretion flow, where the velocity in fig. 2 is now 
a downward velocity, the physical curve is the one with negative 
slope. 
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FIG. 2: Velocity as a function of radius for two different values 
of injection velocity, where the injection surface is at r = 1.2. 
The solid curves show the two relevant solutions to the cubic 
equation (|10p for the case of injection velocity v (1.2) = .34, 
a value very slightly less than the critical value. The heavy 
overlay shows the form of the solution that (when the injection 
is exactly at the critical velocity) connects through the sonic 
point. The dashed curves are for injection at a slightly lower 
velocity, «(1.2) = .3; the lower dashed curve is a physical 
possibility in which the flow remains subsonic in the entire 
space. 



For the moment, we keep the source turned off, again 
taking boundary conditions for injection at a surface, but 
now turning on the injection at t = and solving the 
time dependent equations, to see if there is a stable de- 
velopment of steady flow at sufficiently long times. Our 
computational abilities are not sufficient to deal with the 
sonic transition for the time dependent case, but it is 
illuminating to look at the time behavior both in a com- 
pletely subsonic case and in a completely supersonic case. 

We begin with the case of equations (O, ([6]) with the 
source terms set equal to zero, and consider the two 
generic cases of injection through a surface located at 
some ro > 3/2. Starting at time t = with nearly empty 
space above the injection surface, we maintain a con- 
stant energy density and inflow velocity on the surface. 
To avoid a singular space derivative at the boundary we 
took an initial energy density profile exp— a(r — 1) in- 
stead of absolutely empty space for r > r at t — as 
the initial condition. In figs 2, 3 we show time behavior 
for the supersonic case ro = 7, a — 1.5, and vo = .7. 

We see that a steady flow solution does indeed establish 
itself in short order , in a continually expanding region, 
a flow in which the initial kinetic energy, determined by 
the injection velocity vq, is being converted to internal 
energy. The profiles one sees at the rear of the solid 
figures show steady flow conditions with the velocity plot 
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conforming to the solution of ([TO]) . 




FIG. 3: Energy density, p (or T 4 ), in arbitrary units, for the 
case of steady subsonic injection at a distance ro = 7r s , with 
injection velocity v(l) — .15. The injection is turned on at 
t = 0. The plot shows the rapid evolution into steady flow in 
the region 1 < r/ro < 2, and the conversion of kinetic energy 
to internal energy in the decelerating flow. The flow moves 
from left to right in the figure. 



v 




FIG. 4: Outward velocity for the case of steady subsonic in- 
jection, under the same conditions as used in fig. 3, showing 
the slowing of the flow as it moves upwards. 

For comparison, in figs. 5,6 we show the case of super- 
sonic injection, with initial velocity v — .7. Again, the 
steady flow establishes itself quickly. We see the matter 
speeding up rather than slowing down, with internal en- 
ergy being traded for kinetic energy, as in the standard 
scenario in the literature. 

We have used the same codes to check on the stability 
of solutions against radial perturbations, by taking an 




FIG. 5: Energy density for the case of supersonic injection 
with injection velocity of .7c . 




FIG. 6: Velocity evolution for the case of supersonic injection. 
Taken together with fig. 5, it shows the upward acceleration 
through conversion of internal energy to kinetic energy. 

initial distribution that is a fit to the steady flow limit, 
adding localized perturbations at different distances, and 
confirming rapid relaxation back to the steady flow so- 
lution, both for the cases of subsonic injection and for 
supersonic injection. 

The above discussion was given in order to illuminate 
the problem that we set out to solve, but that prob- 
lem was defined by a given rate of energy-momentum 
deposition in the entire space rather than by injection 
rates through a boundary surface. Turning to the case 
h°,h r 7^ in ([5]), ©, still for the case of steady flow, 
we now ask for solutions for the whole space, from the 
horizon at r = 1 to r = oo. 

We begin with an example that can be reduced to a 
cubic equation by quadrature (although it is both artifi- 
cial and doesn't address the case h r — 0, h° ^ in which 
we are the most interested). This is the case in which 
h r — = 0, so that the right-hand side of © is 

zero, and we can use p — C^ x y~^ to eliminate p in 
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This choice corresponds to injecting pure energy density 
and no momentum density as seen from the frame moving 
with the plasma. Thus it doesn't satisfy the ground rules 
of being externally provided in the observer-at-infinity 
frame, and not dependent on the details of the flow itself. 
But the example gives a neat illustration of the behavior 
of solutions when we have volume injection rather than 
injection through a lower boundary surface. 

Following ref. [H we chose a power-law radial depen- 
dence, 



h°(r) = r- 5 , 
-vf^h = 0, 



(12) 



where in the first of these equations, for the rate of energy 
deposition ho, we use used the exponent suggested by 
ref.[5|. In this reference, the dependence of the source 
strength on position in the region near the black hole 
is taken as (rcos(#))~ 5 (in some limited angular region) 
where 9 is the polar angle. The polar angle dependence 
of the source is important, and we treat some aspects of 
it in the next section, after exploring examples of purely 
radial flow. 

Integrating (|5|) for the source choice lfT2|) we obtain, 



= dCl -1/r) 



(rr 2 -r- 2 ) 



(13) 



where the integration constant r\ is the stagnation ra- 
dius. 

For the plasma to escape at all there must be a sonic 
point r < T\ where, as we move up, we go through the 
value v(ro) = — 1/ All energy injected below the stag- 
nation radius tq goes down the hole. By trial and error, 
we can adjust the two integration constants d and r\ in 
order to find the single solution that passes through the 
two sonic surfaces that now present themselves, that is, 
the transitions at v = ±3 -1 / 2 . In fig. 7 we show basically 
this solution, but with the integration constants detuned 
very slightly so that one sees the normal avoidance of 
crossing of the different solutions of the cubic. 

In this solution approximately one half of the energy 
deposited goes down the hole. There is also a one param- 
eter group of solutions that go through the first sonic 
surface, and pass through a v — surface, but which 
never go supersonic, with the velocities decreasing at 
large heights as described previously. 

In fig. 8 we show a solution for the same equations, but 
with slightly different input parameters d and n which 
lead to a solution where some of the flow does escape, 
but in which the flow remains subsonic and ultimately 
slows. 

As we noted above, to get the second of equations 
(fT2|) we assumed that energy density, and no momen- 
tum, is deposited at given rate in the co-moving system, 
rather than in the observer-at-infinity frame, leading to 
© . This is consistent with the assumptions of ref. [B| , as 




FIG. 7: The solution to (|13[) for the unique parameters that 
connect the supersonic ingoing flow at the horizon to super- 
sonic outgoing flow at large distances. The red curve gives 
the critical path that can be achieved when the parameters 
are exactly tuned. 



0.5 



-0.5 

-1 




FIG. 8: The red curve shows the one parameter set of solu- 
tions to ()13[) in which flow does escape to infinity but remains 
subsonic. The input parameters Ci and ri are only slightly 
different from those used in fig. 7, with the stagnation point 
ri being slightly farther out. 



we can see from reverse engineering eq. 2 of this refer- 
ence, which discusses the relationship of the source term 
in the two frames. 4 In the present case, the assumption 
allows the differential equation for the velocity field to 



4 We note also that h r 
equations used in ref. 



— vf = is required to obtain the 

[l5l in a rather similar application 
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be reduced to a simple cubic algebraic equation, shown 
below. However with this choice the rate of deposition 
of energy-momentum in the plasma depends through v 
on the motion of the plasma, which it should not, in the 
physical case. 

Next we do the parallel analysis for the case of h°(r) — 
r ~ 5 /( r ) h r (r) — 0, that is, pure energy deposition as 
viewed by the observer at infinity. We numerically solve 
the coupled differential equations, in the three regions: 
horizon to first sonic point; first sonic point to second 
sonic point; second sonic point to infinity. Adjusting 
the initial conditions in each segment, by trial and er- 
ror, leads to a plot that is so similar to the heavy curve 
in fig. 7 that we do not need to plot it. The v — point 
is at n = 1.48, and we find that 44% of the energy de- 
posited escapes; the remainder goes down the hole. This 
solution that becomes supersonic is unique, but again 
there is a one-parameter family of perfectly acceptable 
solutions that remain subsonic. This leads us to our first 
important conclusion: Specification of steady energy and 
momentum deposition rates throughout the whole space 
does not determine a unique steady flow solution. 

Presumably, which, if any, of the steady flow solutions 
is obtained a long time after some turn-on of the sources 
will be determined by the solution of a time dependent 
problem, and the solution could depend delicately on pre- 
existing conditions such as the configuration of plasma 
present prior to the turn-on or it could depend on the 
profile of the turning on function. We do not have the 
resources to explore the establishment of the sonic points 
through solution of the time dependent equations, except 
in a fragmentary way. However, in fig. 9 we show the 
results of a time dependent solution where we have both 
a lower bounding surface with fixed subsonic outward 
velocity on the surface, and volume injection of energy 
defined by ft- = r~ 5 ,h r = 0, where the injection is turned 
on at t — 0. In the profile at the rear of the 3D plot we 
see the steady flow solution that develops, which remains 
subsonic. This calculation is supportive of the idea that 
stable, totally subsonic flows are a possible outcome. 

Finally, going back to the steady flow equations, we 
have added massive nucleons, in small enough percent- 
ages so that they contribute significantly to the energy 
density but negligibly to the pressure. Here, in the su- 
personic case, there is an interplay between the injection 
radius and initial nucleon contamination that can limit 
the asymptotic value of T to values much lower than the 
range ~ 100 that is desired. A qualitative statement is 
that when the energy density in the e gas is fifty times 
the rest energy density in the nucleons, the Lorentz fac- 
tor will approach 70 asymptotically, for the case of in- 
jection at 7 times the Schwarzschild radius and at an 
initial speed 5% greater than the sound speed. If we 
take the e ± plasma to have a temperature of 4.5MeV, 
at injection, this corresponds to a nucleon mass density 
Pn ~ 10 6 <7c~ 3 . If we operate with significantly less than 



V 




FIG. 9: An example of a time dependent solution of 
with source term as described in text, above a bounding sur- 
face at r = 1.75, where the injection velocity is .1, showing, 
in the rear profile of the figure, a steady flow that remains 
subsonic over all space above the boundary surface, with the 
velocity reaching a maximum at about r — 2.3 

this amount of contamination, the influence of gravity on 
the flow is negligible. 

The results of this section seem to open the door to a 
modification of the standard picture of the very inner core 
of a v driven jet, in which in the innermost region, where 
energy deposition into the plasma is from the annihilation 
source, the plasma may be decelarating. It is true that 
the authors of ref Q , in a time dependent simulation, find 
a rapid evolution to a steady flow solution that makes 
both sonic transitions without incident. This may be 
completely due to the angular dependence built into the 
authors' z~5 dependence of the source term, or it could 
be because the time-dependent simulation chose a steady 
flow limit embodying the sonic transitions. The plot of T 
shown in fig. 4 of ref. [fj is on such a large scale of r as 
to make the region of neither transition visible in detail. 

Moreover, even though we use a function similar to 
that of these authors for the rate of energy deposition (as 
seen by the observer at oo), it appears that this function 
enters into our equations in a different way. We note 
equation (2) in ref [B| , relating the rate of energy deposit 
q as seen by the observer at oo to the rate q' as seen in 
the comoving system, 

where 7 = 4/3 for our case. This equation correctly de- 
scribes the relation for the case in which the external 



8 



source is taken to have T' r0 = in the comoving sys- 
tem, that is, to be at rest in this system. As mentioned 
above, this makes the rate of momentum deposit by the 
annihilation process dependent on the flow velocity of 
the plasma. Our simulations confirm that the assump- 
tion gives a bias toward escape in the following sense: 
even though we find that specifying the source does not 
uniquely determine the steady flow, if we compare the 
solutions which do pass through both sonic points, and 
accelerate all the way to r — oo, a significantly smaller 
fraction of the deposited energy is captured in the case in 
which the pure energy is added in the comoving system, 
as opposed to the observer-at-infinity system. 

In any case, if we need to go to a time dependent solu- 
tion to pick which of the different steady state flows even- 
tuates (for a fixed energy density injection profile), then 
we would like to know what external factors affect the 
answer. It might be helpful in understanding the close-in 
physics to know what the large relativistic hydrodynam- 
ics codes give for the detailed flow near the hole, the 
location and shape (for non spherical solutions) of both 
sonic surfaces, as a function, say of a turn-on profile for 
the source or the distribution of initial contamination by 
nucleons. 

In any case, in view of our results on the physical plau- 
sibility of the flows that remain subsonic, we assume for 
the remainder of this paper that the (conventional) anni- 
hilation engine indeed produces a plasma that slows (but 
that does not cool much) as it moves outwards, filling the 
close-in regions. This plasma can be helpful in providing 
a seed for the plasma halo that we discussed in the intro- 
duction, but the energetics of this halo will be sustained 
largely by v scattering off of the plasma particles rather 
than by v, v annihilation processes. 

ANISOTROPIC INJECTION 

The considerations of the previous section were con- 
fined to a region extending to a distance of two or three 
times the horizon radius. Now we go to a distance of 
roughly ten times the horizon radius, as sketched in fig. 
1, and discuss the dynamics of a nonisotropic plasma, 
with higher densities at higher angles since the higher 
angle region at a given distance has greater proximity 
to the neutrino-emitting torus. We shall demonstrate a 
path through the sonic transition that depends exactly 
on this angular dependence. Gravity is nearly irrelevent 
at the distances in question. Nor will we include a source 
term, as in (JXJ) , although there will be heating from 
v scattering and cooling, the latter from e + , e~ annihila- 
tion, roughly in balance. 

Here we consider only the steady flow case. But the 
time dependent solutions presented for the isotropic case 
give us confidence that steady flow solutions will estab- 
lish themselves in a growing volume above the injection 



point as time goes on. We begin again with JT]) , but now 
take p(ro, 9) = const. (I + .7 smO) on the inner boundary. 
This corresponds to a temperature that is roughly 10% 
greater at the equator than at the pole. In figs. 8,9 we 
show an example of the resulting flow as determined from 
a solution in a region 6 < 9 < it/2, r < r < r\. The 
solution in the region tt/2 < 9 < tt, would be obtained 
by up-down symmetry. The derivatives 8/ 89 of all quan- 
tities vanish for 9 = tt/2. In the example shown we took, 
Oo = -12 



1 . 2 
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FIG. f0: Density contours in the steady flow solution for 
the non-isotropic case, showing the increases of magnitude 
in both components of the velocity as we approach the sonic 
singularity of the near corner of the plot. This point is just 
one point on the sonic surface that we explore by repeatedly 
using different values of (ri,0i) (the corner coordinate), each 
of the pairs being chosen to put the maximum velocity point 
on the corner. 




FIG. f f : Radial and tangential velocity contours in the steady 
flow solution for the non-isotropic case. The upper sheet 
shows the radial velocity, and the lower the tangential ve- 
locity, with negative values since it is inwards. At the near 
corner of the plot, ri, 6\ 

We designate the nearest corner on the lower plane of 
the bounding cubes in figs. 9,f0 as the the point r\,9\. 
Approaching this corner, we see a rapid rise in v r as well 
as the increase in the magnitude of v t = (which is 
negative). In the example shown, the resultant speed at 
this corner is very close to the speed of sound (l/y3)- 
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Repeating this calculation with different input values of 
9q allows us to plot out the envelope of the supersonic 
region. 

Going back to fig.l we show the part of the sonic 
boundary that is nearest to the injection surface as a pair 
of curves sprouting out of the top of the injection surface. 
As we noted, the velocity field has a gradient approaching 
infinity as we approach the surface from the outside, with 
the derivative of the energy density approaching (nega- 
tive) infinity as well, as it must to conserve energy. On 
the inside of the surface, as we we move upwards we then 
will have the continuation of upward (but now decreas- 
ing) acceleration to velocities approaching the speed of 
light. 

In the previous section we investigated the sonic tran- 
sition in isotropic flows. These transitions were driven 
by a combination of gravitational effects and the radial 
dependence of the source terms; in the absence of both 
there would have been no sonic transition. 

In the present section we found that a non-isotropic 
flow can make the transition on its own, that is to say, in 
the absence of gravitational effects or infusion of energy 
from the outside. This encourages us to take the model 
described in the introduction seriously. In this model 
the plasma halo moves outward with small velocities in 
a region in which v, v annihilation no longer provides 
a big energy source, the heating being instead through 
v, e ± collisions. The anisotropic temperature created by 
this heating mechanism gives the needed precondition for 
a more or less conical sonic surface to develop, and for 
accelerating supersonic flow to prevail in the interior of 
the cone, was shown schematically in fig. 1. 5 



DISCUSSION 

We search for a mechanism that can accelerate plasma 
to supersonic speeds, and thence to a Lorentz factor of 
up to r = 100. The second stage appears to be easier 
than the first; our solutions for a spherically symmetri- 
cal case with isotropic supersonic injection are similar in 
radial structure to the models of refs [|| [Bj], where the 
injection is confined to a cone (and lateral spreading is 
restricted by pressure from the medium.) Even in these 
models with a confined jet it should be the case that the 
plasma configuration is more spherical than jet-like in 
the innermost region, r < (a few)r s . Furthermore, once 
we enter the domain v ~ c we recover the scaling rela- 
tion r(r) ~ r found in the analytic results of ref |3j|, and 
confirmed approximately in the numerical results of ref. 



5 The authors of ref. [IH , after analyzing energy deposition by neu- 
trinos, found higher temperatures off-axis than along the axis and 
suggested that this might strengthen the jet 



5]. However, in the spherically symmetrical case with 
subsonic injection, the flow slows as it moves upwards, 
stagnating at r = oo. 

The situation changes when we introduce a dependence 
of the injection energy density on the polar angle 9. Dis- 
tinguishing two cases, roughly, we can consider either 
an injection temperature that is greater at the equator, 
the case that corresponds to our qualitative picture of 
how the halo forms, or an injection temperature that is 
greater at the poles. In the first case we indeed find a 
flow that accelerates in the region near the axis, with the 
resultant of radial and tangential velocities reaching the 
sound speed as we approach, from the outside, a hornlike 
surface (that is, cone-like with expanding angle), centered 
on the axis. 

The hydrodynamical equations are singular as we ap- 
proach this surface. For the boundary condition for the 
region on the inner side of this surface we take the values 
of energy density and directions of velocity from the pre- 
vious calculation, boosting the magnitude of the velocity 
to a value barely above the speed of sound, providing a 
new boundary condition for the acceleration up to large 
r in the inside of the sonic cone. 

In the case in which the subsonic injection density is 
greater toward the equator than toward the poles, we find 
that the velocity decreases, going upward, for all values 
of 9. We conclude that having the injection temperature 
higher toward the equator is essential in forming the jet 
in the vicinity of the polar axis. 

All of our simulations in the non-isotropic case were for 
a completely relativistic plasma, P = p/3. But we be- 
lieve that the numerical experiments that we mentioned 
in sec. 2 for a plasma with nucleoli loading support the 
conclusion that in the domains of density of p < 10 6 <?c -3 
the qualitative behavior is not changed by the loading. 

It would be interesting to see the detailed results of 
large time dependent codes that have been applied to 
the study of the large scale structure of the jets, for ex- 
ample, GENESIS, applied locally in the region near the 
black hole under a number of assumptions for the initial 
conditions that prevail, in order to better understand the 
universe of possibilities in this region 

Comparing to the model of ref. which assumed a 
disk with a sharp boundary in which the density is inde- 
pendent of height above the median plane, we find that 
even if we begin with this system, and a (near) vacuum 
on top, a pair-plasma will be created in the region above 
the ^-surface, at a temperature somewhat less than that 
of the ^-surface. The seed would be from D + v — > e + + e~ 
interactions, either produced locally or drifted in from be- 
low, but the heating to final temperatures is dominantly 
from v — e ± collisions. Estimates leading to these con- 
clusions are given in the appendix. 

This hypothetical time-dependent scenario deserves to 
be addressed with a real transport calculation, which we 
have not done. However, we had previously looked [l2T | 
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at a steady flow solution for the the vertical structure of 
an atmosphere with a sharp surface, with accretion rates 
such that the optical thickness (to z/s) is large, where nu- 
clear opacity, relativistic opacity, and vertical hydrostatic 
equilibrium were taken into account. We found that the 
nuclear densities must decrease steadily as we go up, and 
end at a very small value at the surface, a finding that 
coexists well with the assumptions of the present paper. 

That said, we return to reasons why the present work 
is not in detail a realistic model of the jet formation: 

1) We may have been over-optimistic in taking small 
nucleon densities (p < 10 8 ) in the region of larger angles 
9. But this is the problem that we could formulate and 
solve. Of course, as we approach 9 = n/2 the densities 
get really large. However, had we constrained to a cone 
bounded by 9 max we would have had to put an arbitrary 
boundary condition on the surface of the cone, whereas 
in the solution that has the full 4ir solid angle, there 
is no such boundary, only the up-down symmetry that 
guarantees that the derivatives dj 89 vanish at 9 — ir/2. 

2) Our injection sphere had an arbitrary radius, so the 
curve which in fig. 1 shows the sonic surface is thus ar- 
bitrarily situated as well. (Of course, our assumption 
that the injection velocity is radial is arbitary as well.) 
Suppose that we had taken the radius to be a little bit 
smaller, and then calculated the temperatures and veloc- 
ities on the first surface, as a new boundary condition. 
Then transverse velocities would no longer vanish on the 
injection surface. We would then have found the horn- 
shaped sonic surface to start lower, in fact lower than 
the injection radius for a 6 = 0. Our actual choice of 
7r s (sa 50km, for a three solar mass BH) was motivated 
by a compromise between these considerations and the 
need both to get outside of regions of high nucleon den- 
sity and regions in which the back reaction from neutrino 
transport would be a large perturbation on the perfect 
fluid dynamics. 

For the moment ignoring these drawbacks, we ask how 
the energy injected in the outward beam by the mech- 
anism of the present paper could compare with that in, 
e.g., the models of ref.[5j. In a representative example of 
the latter models, as the jet passes our surface at 50 km., 
it carries something like 10 50 — 10 51 ergs/sec. across this 
boundary. The temperature of the plasma in the jet is 
already well less than 1 MeV at this point. By contrast, 
in our guesswork example, the temperature just above 
the injection radius, but at somewhat higher angle 9 is 
of the order 3 MeV. As is clear from fig. 6, as we move 
to smaller values of theta and approach the sonic sur- 
face, the energy density (or T 4 ) drops rapidly, and it will 
continue to drop as we move away from the sonic sur- 
face in the inside supersonic region. However the plasma 
energy density inside the lower part of the sonic "horn" 
still should be greater in this model than in the model of 
ref. @ 

In this more muscular jet, the tolerance for baryonic 



contamination would likewise be much greater. From the 
standpoint of total energetics, the new scenario uses the 
ability of the plasma halo, which has greater volume than 
the farther- in annihilation region of ref [3. to capture 
energy from a larger fraction of the total v output of 
the torus and deliver some fraction to the jet. Some 
more detailed considerations on the energetics of heating 
mechanisms are given in the appendix. 

The core results of this paper remain the solutions of 
the hydrodynamical equations in the perfect fluid do- 
main. Some of the qualitative properties of these so- 
lutions have not surfaced in other results of numerical 
computations known to us. At the least, it would be in- 
teresting to test the more comprehensive codes against 
our results, or against similar computational fragments 
based on simplified systems. As we learned during the 
generations of supernova calculations, once "all the im- 
portant physics" is put into calculations of two different 
groups, it is certain that the calculations are not exactly 
the same as each other, and comparison of the results is 
problematical. 
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APPENDIX 

Here we discuss the energetics of the configuration 
sketched in fig. 1, which provided some of the motiva- 
tion for the considerations of sec. 3. The figure shows a 
slice of the system by a plane that is perpendicular to the 
plane of the accretion disk and which contains the axis of 
the disk. The cross-sections of an inner torus and that of 
an outer torus are shown with different shadings. In the 
inner one we expect the neutrino opacity to be dominated 
by scattering and absorption on nucleons, as it is in the 
simulations of refs [1] The nucleon density decreases 
steadily as one moves from the center to the edge of this 
torus; at the edge completely leptonic opacity is becom- 
ing important. In the surrounding toroidal region, less 
shaded, as we move from the inside to the outside the 
leptonic opacity becomes dominant. For our choice of 
parameters the v free path in this region is about 30 km. , 
a little too long to signal v trapping, but short enough 
to be relevant. For a case with slightly higher tempera- 
ture we would have actual trapping. Thus we picture the 
^-surface as lying somewhere between the inside and out- 
side of the region with medium shading. Outside of this 
region, in the very lightly shaded region, the v opacities 
are small enough so as not to change the v fluxes and 
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energy distributions at the surface of the circumscribed 
sphere, our "injection" sphere. 

A shell of plasma at the edge of the lightly shaded re- 
gion, for example at the point B labeled "T = 2.9 MeV" 
is cooled by thermal v emission and heated by scattering 
of neutrinos from below. The latter have an energy dis- 
tribution corresponding a ^-surface at T = 5 MeV, but 
their number density is reduced at the point B by a fac- 
tor X~(solid angle subtended at B by v - surface) / Ait w .1. 
The temperature at point B will be set by balancing the 
heating and cooling rates. 

To estimate the rate of cooling we begin with the total 
cross-section for the process e~ + e + — > v e + D e , in the 
completely relativistic limit, with the four- momenta of 
the electrons denoted by p, q, and taking sin 2 (9w) = -25, 



c e +.e- 



5G|n 
12tt 



(15) 



Defining, Tp as the temperature of the plasma at point 
B, and Tg as the temperature of the neutrino-surface, we 
write the cooling rate, per unit volume, at B as, 



Qcool 



4(2^)" 6 / d 3 pd 3 qf(E p , T B )f(E q , T B ){E p + E q ) 



x (j e + e _ — *3J5Gj?Tj2 



re 3.8 x 10 33 T 1 9 1 ergs (cm)" 3 sec" 1 , (16) 

where f(E,T) is the Fermi distribution with zero chemical 
potential. 

Next we estimate the net heating rate at B beginning 
from the sum of elastic cross-sections, 

°T = ( a u it e- + a Pi,e- + a u it e+ + °P 4 ,e+) 



UG 2 F 
3tt 



-P/*<r 



(17) 



where the sum is over the three v flavors, and where p 
and q are now the respective momenta of the incident 
neutrino and electron. We estimate the rate of energy 
delivery as, 



g hcat = A(2vr)- 6 J d 3 pd 3 q(E p - E q )f(E p , T B ) 
xf(E g ,T s )a T = 3AXG%T^, 



(18) 



where the initial factor of 2 is from the sum over 
spins (the cross-sections quoted are averaged over initial 
spins. The expression (|18p vanishes when T$ = T B , as 
it should. To obtain the last equality in (fl"8"|) we took 
T s = 2^ T B . For this choice, T s = 5MeV T B = 2.9, 
with A rs .1, comparing with (| 16[> . we see that the heating 
from scattering is roughly in balance with cooling from 
e + ,e~ annihilation. 



We check our assertion that the heating through the 
scattering process is greater than that from the conven- 
tional neutrino annihilation process. The rate for the 
latter is, from lfTSl) . 



<?a„n = A 2 4772a e+ , e -(27r)- 6 J d 3 pd 3 qf(E p ,T s )f(E q ,T s ) 

x (E p + E q )a e+ , e - = b.7\ 2 r,G 2 F T% , (19) 

where r\ is a colinearity reduction factor rj — 
(p fJ, q )Ji /(E p E q )), which is about .2 on the axis, at the injec- 
tion radius, in the sketch of fig. 7, and is rapidly decreas- 
ing as 9 increases. Here the pref actor 4 comes from the 
relation of cross-sections, o v+9 ^ e + +e - = 4<r e +_|_ e -_ >l/+p . 
Taking A = .1, we find the rate is much less than the 
energy deposition rate in the scattering reactions (fTSf . 

Finally we estimate the rate of radial momentum trans- 
fer to the outer layer through the scattering process as 
w 3x the rate of energy transfer. With the parameters 
at hand this leads to rough balance with the oppositely 
directed rate of momentum deposition due to gravity 

In this paper, we have more than once remarked on the 
role of z/-lepton reactions as a source of neutrino opacity 
that can affect the location and temperature of the v 
surface. Since this issue appears not to have a qualitative 
effect on our results, we have not developed it in detail. 
Here we note that ref. [4j and others include only the 
lcptonic reaction included in the opacity was v e + v e — > 
e + + e~. However, in conditions in which all three favors 
of v, v are trapped one needs to take into account all 
(zz+lepton) reactions, and for a v of any flavor there are 
11 separate reactions. As it turns out, leptonic trapping 
does not quite occur in the present scenario although it 
could in extreme cases of energy output of the disk. 

We have also looked at the same set of processes as a 
means of establishing the plasma halo in the first place 
(or to put it another way, as a means of heating a re- 
gion in which p nuc < 10 8 , so that reactions on nucleons 
are inconsequential). We find that through the chain of 
first creating some pairs from the annihilation process 
after which the dominant e — v heating takes over, tem- 
peratures approaching the steady state temperature are 
reached promptly. 
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